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Abstract. The evolution of the particle distribution func- 
tions inside a relativistic jet containing an e~-e + pair plasma 
and of the resulting 7-ray and X-ray spectra is investigated. 
The first phase of this evolution is governed by heavy ra- 
diative energy losses. For this phase, approximative expres- 
sions for the energy-loss rates due to inverse-Compton scat- 
tering, using the full Klein-Nishina cross section, are given as 
one-dimensional integrals for both cooling by inverse-Compton 
scattering of synchrotron photons (SSC) and of accretion disk 
photons (EIC). 

We calculate instantaneous and time-integrated 7-ray spec- 
tra emitted by such a jet for various sets of parameters, de- 
ducing some general implications on the observable broadband 
radiation. Finally, we present model fits to the broadband spec- 
trum of Mrk 421. We find that the most plausible way to ex- 
plain both the quiescent and the flaring state of Mrk 421 con- 
sists of a model where EIC and SSC dominate the observed 
spectrum in different frequency bands. In our model the flar- 
ing state is primarily related to an increase of the maximum 
Lorentz factor of the injected pairs. 
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1. Introduction 

Accretion of matter onto a central black hole is the most rele- 
vant process to power active galactic nuclei (Lynden-Bell 1969, 
Salpeter 1969, Rees 1984). However, the details of the conver- 
sion processes of gravitational energy into observable electro- 
magnetic radiation are still largely unknown. The discovery of 
many blazar-type AGNs (Hartman et al. 1992, Fichtel et al. 
1993) as sources of high-energy gamma-ray radiation dominat- 
ing the apparent luminosity, has revealed that the formation of 
relativistic jets and the acceleration of energetic charged parti- 
cles, which generate nonthermal radiation, are key processes to 
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understand the energy conversion process. Emission from rela- 
tivistically moving sources is required to overcome gamma-ray 
transparency problems implied by the measured large lumi- 
nosities and short time variabilities (for review see Dermer & 
Gehrels 1995). 

Repeated gamma-ray observations of AGN sources have 
indicated a typical duty cycle of gamma-ray hard blazars of 
about 5 percent, supporting a "2-phase" model for the central 
regions of AGNs (Achatz et al. 1990, Schlickeiser & Achatz 
1992). According to the 2-phase model the central power- 
house of AGNs undergoes two repeating phases: in a "quiescent 
phase" over most of the time (~95 percent) relativistic charged 
particles are efficiently accelerated in the central plasma near 
the black hole, whereas in a short and violent "flaring phase" 
the accelerated particles are ejected in the form of plasma blobs 
along an existing jet structure. 

We consider the acceleration of charged particles during 
the quiescent phase. The central object accretes the surround- 
ing matter. Associated with the accretion flow is low- frequency 
magnetohydrodynamic turbulence which is generated by vari- 
ous processes as e.g.: 

(a) turbulence generated by the rotating accretion disk at large 
eddies and cascading to smaller scales (Galeev et al. 1979); 

(b) stellar winds from solar-type stars in the central star cluster 
deliver plasma waves to the accretion flow; 

(c) infalling neutral accretion matter becomes ionized by the 
ultraviolet and soft X-ray radiation of the disk. These pick-up 
ions in the accretion flow generate plasma waves by virtue of 
their streaming (Lee & Ip 1987); 

(d) if standing shocks form in the neighbourhood of the cen- 
tral object they amplify any incoming upstream turbulence in 
the downstream accretion shock magnetosheath (McKenzie & 
Westphal 1969, Campeanu & Schlickeiser 1992). 

These low-frequency MHD plasma waves from the accre- 
tion flow are the source of free energy and lead to stochastic 
acceleration of charged particles out of the thermal accretion 
plasma. 

The dynamics of energetic charged particles (cosmic rays) 
in cosmic plasmas is determined by their mutual interaction 
and interactions with ambient electromagnetic, photon and 
matter fields. Among these by far quickest is the particle- wave 
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interaction with electromagnetic fields, which very often can 
be separated into a leading field structure F and superposed 
fluctuating fields SF. Theoretical descriptions of the trans- 
port and acceleration of cosmic rays in cosmic plasmas are 
usually based on transport equations which are derived from 
the Boltzmann-Vlasov equation into which the electromagnetic 
fields of the medium enter by the Lorentz force term. The 
quasilinear approach to wave-particle interaction is a second- 
order perturbation approach in the ratio ql = (5F/F ) 2 and 
requires smallness of this ratio with respect to unity. In most 
cosmic plasmas this is well satisfied as has been established 
either by direct in-situ electromagnetic turbulence measure- 
ments in interplanetary plasmas, or by saturation effects in the 
growth of fluctuating fields. Nonlinear wave-wave interaction 
rates and/or nonlinear Landau damping set in only at appre- 
ciable levels of (5F) 2 and thus limit the value of <?l < 1. We 
assume the AGN plasma to have very high conductivity so that 
any large-scale steady electric fields are absent. We then con- 
sider the behaviour of energetic charged particles in a uniform 
magnetic field with superposed small-amplitude (SB) 2 <C B 2 
plasma turbulence (5E, SB) by calculating the quasilinear cos- 
mic ray particle acceleration rates and transport parameters. 
This is by no means trivial since especially for the interac- 
tion of non-relativistic charged particles with ion- and electron- 
cyclotron waves thermal resonance broadening effects are par- 
ticularly important (Schlickeiser & Achatz 1993, Schlickeiser 
1994). The acceleration rates and spatial transport parameters 
are then used in the kinetic diffusion-convection equation for 
the isotropic part of the phase space density of charged par- 
ticles F(x,p,t) which for non-relativistic bulk speed u <§C c 
reads 

at b °-dx lK dx l 

dF pdudF 1 d , 2 A dF 2 . „ 

Here x denotes the spatial coordinate along the ordered 
magnetic field, p the cosmic ray particle momentum, k is the 
spatial diffusion coefficient, A the momentum diffusion coeffi- 
cient, and So denotes the "Stossterm" describing the mutual 
interaction of the charged particles and their injection. 

With respect to the generation of energetic charged par- 
ticles, the basic transport equation (1) shows that stochas- 
tic acceleration of particles, characterized by the accelera- 
tion time scale tA = p 2 /A, competes with continuous energy 
loss processes pi oss , characterized by energy loss time scales 
th = p/\pi oss \- Dermer et al. (1996) have recently inspected the 
acceleration of energetic electrons and protons in the central 
AGN plasma by comparing the time scales for stochastic ac- 
celeration with the relevant energy loss time scales. At small 
proton momenta the Coulomb loss time scale is extremely sen- 
sitive to the background plasma density and temperature, and 
for slight changes in the values of these parameters cosmic ray 
protons may not be accelerated above the Coulomb barrier. 
Although at small particle momenta the plasma wave's dissi- 
pation and the interaction with the cyclotron waves become 
decisive and might modify the acceleration time significantly, 
the results of Dermer et al. (1996) demonstrate that reasonable 
central AGN plasma parameter values are possible where the 
low-frequency turbulence energizes protons to TeV and PeV 
energies where photo-pair and photo-pion production are effec- 
tive in halting the acceleration (Sikora et al. 1987, Mannheim 



& Biermann 1992). According to the results of Dermer et al. 
(1996) it takes about ~ lOJWs days for the protons to reach 
these energies, where Ms is the mass of the central black hole 
in units of 10 8 Mq. The corresponding analysis for cosmic ray 
electrons shows that the external compactness provided by the 
accretion disk photons (Becker & Kafatos 1995) leads to heavy 
inverse Compton losses which suppress the acceleration of low- 
energy electrons beyond Lorentz factors of 7 « 10 — 100. It 
seems that due to their much smaller radiation loss rate cos- 
mic ray protons are effectively accelerated during the quiescent 
phase in contrast to low energy electrons. 

Now an important point has to be emphasized: once the 
accelerated protons reach the thresholds for photo-pair ("f P: th — 
m e c 2 1 < e >— 5 • 10 4 e 1 _1 ) and photo-pion production and the 
threshold for pion production in inelastic proton-matter col- 
lisions they will generate plenty of secondary electrons and 
positrons of ultrahigh energy which are now injected at high 
energies (E s > 25GeVe^ 1 ) into this acceleration scheme. 
< e >— 10ei eV denotes the mean accretion disk photon en- 
ergy. It is now of considerable interest to follow the evolution 
of these injected secondary particles. 

Although many details of this evolution are poorly under- 
stood, it is evident that the further fate of the secondary par- 
ticles depends strongly on whether they find themselves in a 
compact environment set up by the external accretion disk, or 
not. As has been pointed out by Dermer & Schlickeiser (1993b) 
as well as Becker & Kafatos (1995) the size of the gamma-ray 
photosphere (where the compactness is greater unity so that 
any produced gamma-ray photon is pair-absorbed) is strongly 
photon energy dependent. The gamma-ray photosphere attains 
its largest size at photon energies E v ~ SOe^ 1 GeV. Secondary 
particles within the photosphere having energies E s < E p will 
initiate a rapid electromagnetic cascade which has been stud- 
ied by e.g. Mastichiadis & Kirk (1995), which might even lead 
to runaway pair production and associated strong X-ray flares 
(Kirk & Mastichiadis 1992), and/or due to the violent effect of 
a pair catastrophy (Henri & Pelletier 1993) ultimately lead to 
an explosive event and the emergence of a relativistically mov- 
ing component filled with energetic electron-positron pairs. 

In contrast, if the secondary particles are generated out- 
side the gamma-ray photosphere the secondary electrons and 
positrons will quickly cool by the strong inverse Compton losses 
generating plenty of gamma-ray emission up to TeV energies. 
As solutions of the electron and positron transport equation (1) 
for this case demonstrate (Schlickeiser 1984, Pohl et al. 1992) 
a cooling particle distribution (N(p) oc p~ 2 ) with a strong cut- 
off at low (but still relativistic) momentum p c develops, which 
grows with time as more and more protons hit the photo-pair 
and photo-pion thresholds. Such bump-on-tail particle distri- 
bution functions, which are inverted (dF/dp > 0) below p c , 
are collectively unstable with respect to the excitation of elec- 
tromagnetic and electrostatic waves such as oblique longitudi- 
nal Langmuir waves (Lesch et al. 1989). As described in de- 
tail by Lesch & Schlickeiser (1987) and Achatz et al. (1990), 
depending on the local plasma parameters (mainly the elec- 
tron temperature of the background gas and the density ra- 
tio N /n e of relativistic electrons and positrons to thermal 
electrons) these Langmuir waves either (1) lead to quasilin- 
ear plateauing of the inverted distribution function and rapid 
collective thermalization of the electrons and positrons, or (2) 
are first damped by nonlinear Landau damping, but ultimately 
heat the background plasma strongly via the modulation in- 



3 



stability once a critical energy density in Langmuir waves 
W c = n e (kBT e /m e c 2 ) 2 has been built up. Both relaxation 
mechanisms terminate the quiescent phase of the acceleration 
process. The almost instantaneous increase of the background 
gas entropy due to the rapid modulation instability heating 
again leads to an explosive outward motion of the plasma blob 
carrying the relativistic particles away from the central object. 

As we have discussed, in both cases it is very likely that 
at the end of the quiescent phase an explosive event occurs 
that gives rise to the emergence of a new relativistically mov- 
ing component filled with energetic electron-positron pairs. It 
marks the start of the flaring phase in gamma-ray blazars. The 
initial starting height of the emerging blob z% entering the cal- 
culation of the gamma-ray flux should be closely related to the 
size of the acceleration volume in the quiescent phase mainly 
determined by the maximum size of the gamma-ray photo- 
sphere. This scenario is supported by the measurements of 
Babadzanhanyants &: Bclokon (1985) that in 3C 345 and other 
quasars optical bursts are in close time correlation with the 
generation of compact radio jets. A similar behaviour has also 
been observed during the recent simultaneous multiwavelength 
campaign on 3C 279 (Hartman et al. 1996). Further corrobo- 
rative evidence for this scenario is provided by the recent dis- 
covery of superluminal motion components in the gamma-ray 
blazars PKS 0528+134 (Pohl et al. 1995) and PKS 1633+382 
(Barthel et al. 1995) that demonstrate a close physical con- 
nection between gamma-ray flaring and the ejection of new 
superluminal jet components in blazars. 

It is the purpose of the present investigation to follow the 
time evolution of the relativistic electrons and positrons as the 
emerging relativistic blob moves out. Because of the very short 
radiative energy loss time scales of the radiating electrons and 
positrons it is important to treat the spectral evolution of the 
radiating particles self-consistently. In earlier work Dermer & 
Schlickeiser (1993a) and Dermer et al. (1997) have studied the 
spectral time evolution from a modified Kardashev (1962) ap- 
proach by injecting instantaneously a power law electron and 
positron energy spectrum at height Zi at the beginning of the 
flare, and calculating its modification with height in the rel- 
ativistically outflowing blob due to the operation of various 
continous energy loss processes as inverse Compton scatter- 
ing, synchrotron radiation, nonthermal bremsstrahlung emis- 
sion and Coulomb energy losses. Here we generalize their ap- 
proach by accounting for Klcin-Nishina effects and including as 
well external inverse-Compton scattering as synchrotron self- 
Compton scattering self-consistently. 

The acceleration scenario described above leads us to the 
following assumptions on the distribution of pairs inside a new 
jet component at the time of their injection: If the pairs are 
created by photo-pair production, their minimum Lorentz fac- 
tor is expected to be in the range of the threshold value of the 
protons' Lorentz factor for photo-pair production. We use the 
standard accretion disk model by Shakura & Sunyaev (1973), 
which we will describe in more detail in the next section, to fix 
this threshold, determined by the average disk photon energy 
(e). The pair distribution above this cutoff basically reflects the 
acceleration spectrum of the protons, i. e. a power-law distribu- 
tion with spectral index 2<s<3 (71(7) ~ 7~ s ) which extends up 
to 72± ~ 10 6 . This yields the initial pair distribution functions 



with 4<a<5. The differential number of particles in the 
energy intervall [7±,7± + chf±] per unit volume is then given 
by 



d n± = 4 7T N± 7 ± /3± f± (7±)<*7± 



(3) 



where N± is a normalization factor related to the total 
particle density n± through N± « n ± 2 ^'7i±' i ■ 

The detection of TeV 7-rays from Mrk 421 suggests that 
such components must be produced/accelerated outisde the 7- 
ray photosphere for photons of energy ~ 1 TeV. The height of 
this photosphere due to the interaction of 7-rays with accretion 
disk radiation will be determined self-consistently in section 
4. Backscattering of accretion disk radiation by surrounding 
clouds is negligible in the case of BL Lac objects emitting TeV 
7-rays (Bottcher & Dermer, 1995). 

We point out that most of our basic conclusions are also 
valid if the pairs inside the blob are accelerated by other mech- 
anisms, e. g. by a relativistic shock propagating through the 
jet. 

Beginning at the injection height (henceforth denoted as 
Zi), we follow the further evolution of the pair distribution and 
calculate the emerging photon spectra. 




Fig 



1. Model for the geometry of a relativistic AGN jet 



/±(7±)=7i° 



7i± < 7± < 72± 



(2) 



The negligibility of pair absorption due to the interaction 
with the synchrotron and 7-ray emission from the jet is checked 
self-consistently during our calculations. 

Interactions of the jet pair plasma with dilute surrounding 
material will cause turbulent Alfven and Whistler waves. It has 
been shown by Achatz & Schlickeiser (1993) that a low-density, 
relativistic pair jet is rapidly disrupted as a consequence of the 
excitation of such waves. Thus, to insure stability of the beam 
over a sufficient length scale, we need that the density of pairs 
inside the jet exceeds the density of the surrounding material. 
In this case, pitch angle scattering on plasma wave turbulences 
leads to an efficient isotropization of the momenta of the pairs 
in the jet without destroying the jet structure. Thus, additional 
assumptions on our initial conditions are that the particle mo- 
menta are isotropically distributed in the rest frame of a new 
jet component (blob) and that the density of pairs in the jet 
rij 3> n b where n b is the density of the background material. 
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This study is devided into two papers. In the first (present) 
paper we investigate the details of electron/positron cooling 
due to inverse-Compton scattering and follow the pair distri- 
bution and photon spectra evolution during the first phase in 
which the system is dominated by heavy radiative losses. In the 
second paper we will consider the later phase of the evolution 
where collisional effects (possibly leading to thermalization) 
and reacceleration become important, and a plausible model 
for MeV blazars which follows directly from our treatment will 
be presented (Bottcher, Pohl & Schlickeiser, in prep.). 

In section 2 of this first paper, we describe in detail how 
to calculate the energy-loss rates due to the various processes 
which we take into account and give useful approximative ex- 
pressions for the inverse-Compton losses (as well scattering of 
accretion disk radiation as of synchrotron radiation) , including 
all Klein-Nishina effects. In section 3, we discuss the relative 
importance of the various processes. The location of the 7-ray 
photosphere for TeV 7-rays is briefly outlined in section 4. The 
technique used to follow the evolution of the pair distributions 
is described in section 5. In section 6, we describe how to use 
the pair distributions resulting from our simulations in order 
to calculate the emanating 7-ray spectra, and in section 7 we 
discuss general results of our simulations giving a prediction 
for GeV - TeV emission from 7-ray blazars which due to the 
lack of sensitivity of present-day instruments in this energy 
range could not be observed until now. Only two extragalactic 
objects have been detected as sources of TeV emission, namely 
Mrk 421 (Punch et al. 1992) and Mrk 501 (Quinn et al. 1995). 
In Section 8, we use our code to fit the observational results 
on the broadband emission during the TeV flare of Mrk 421 in 
May 1994 and on its quiescent flux. We summarize in section 
9. 



2. Energy- loss rates 

We first consider in detail the various processes through which 
the pairs inside a new AGN jet component lose energy. The 
blob is assumed to move outward perpendicularly to the ac- 
cretion disk plane with velocity c /3r = c y/l — 1/T 2 where T 
is the Lorentz factor of the bulk motion. The mechanisms that 
we take into account are inverse-Compton scattering of accre- 
tion disk photons, synchrotron and synchrotron-self-Compton 
(SSC) losses. It is well-known that energy exchange/loss due to 
elastic (M0ller and Bhabha scattering) and inelastic scattering 
(pair bremsstrahlung emission) do not contribute significantly 
for ultrarelativistic particles. We consider them in detail in the 
second paper of this series, dealing with mildly-relativistic pair 
plasmas. The same is true for pair annihilation losses. 
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In general, the energy loss rate due to inverse-Compton 
scattering is given by the manifold integral 
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Fig. 2. Definition of the angles for calculation of the in- 
verse-Compton losses. 'L' in the right panel denotes the motion of 
the labor frame with respect to the electron rest frame. The subscript 
'ph' refers to quantities of the photon 



fi — cos ip = rjph r] e + sin 6^ sin 6 e cos <f> 
fi' s = cos ip' s = k ji + sin x sin V*' cos 0s • 
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(7) 



Here, primed quantities refer to the rest frame of the elec- 
tron, and the subscript 's' denotes quantities of the scattered 
photon. The definition of the angles is illustrated in Fig. 2. 
Now, let the superscript V denote quantities measured in the 
rest frame of the accretion disk. Then, the differential number 
of accretion disk photons in the blob frame is 
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2.1. External inverse-Compton losses 

We are now considering the single-particle energy loss rate due 
to inverse-Compton scattering of external photons coming di- 
rectly from a central source, which we assume to be an ac- 
cretion disk. We use the accretion disk model of Shakura & 
Sunyaev (1973) predicting, for a central black hole of M ~ 10 6 
- 10 8 Mq , a blackbody spectrum according to a temperature 
distribution T(R) given by 
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Using this photon number and the full Klein-Nishina cross 
section, Eq. (5) becomes 
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where 

F = 1 + £7(1-^)(1-k'). 
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and i?min and i?max are the radius of the inner and outer 
edge of the accretion disk, respectively. 

If all scattering occurs in the Thomson regime (F = 1), we 
find (neglecting terms of order I/7 2 ): 
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where x = \/R 2 + z 2 . Fig. 3 shows the energy loss rate 
computed using the full Klein-Nishina cross-section compared 
to the calculation in the Thomson limit as quoted above as 
well as the calculation in the Thomson limit, combined with a 
point-source approximation for the accretion disk (e. g. Dermer 
& Schlickeiser 1993). 
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Fig. 3. Energy loss rates of a test electron/positron due to in- 
verse-Compton scattering of accretion disk photons. Total disk lu- 
minosity: Ln = 10 44 erg s _1 , height: z = 2 ■ 10 — 4 pc 

For electron energies of 7 <JC Y/(e diak ), Eq. (12) is a very 
good approach to the Inverse-Compton losses. In the case of 
small distances to the accretion disk (z<\0~ 2 pc in the case of 
a disk luminosity L = 10 46 erg s" 1 ; z<10~ 3 pc for L = 10 44 
erg s _1 ) the point-source approximation is not an appropriate 
choice. 



A very useful approximation for all electron energies is 
based on replacing the integration ± J dr] e by setting r) e = 
0. Furthermore, one can approximate the thermal spectrum, 
emitted by each radius of the disk, by a 8 function in energy, 
n(e, Q) « AS(e- < e > [R]) where 
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With these simplifications, and neglecting terms of order 
I/7 2 < 1, Eq. (9) becomes 
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For small values of E, one should use the Taylor expansion of 
Eq. (15), namely 



= 2tt 




-f E (l + la 2 )+f E 2 (l + 3a 2 +a 4 )j+o(V). (16) 

This result (using Eq. [16] for E < 0.05) is illustrated by 
the dot-dashed curve in Fig. 3. In Fig. 5, the inverse-Compton 
losses (according to Eq. [13]) for set of parameters which is 
assumed to be typical for BL Lac objects (Lq — 10 44 ergs _1 , 
2 = 2- 10~ 4 pc, r = 10) are compared to the synchrotron and 
SSC energy losses derived in the following subsections. 



2.2. Synchrotron and synchrotron- s elf- Compton losses 

If, for the synchrotron (sy) losses, we neglect inhomogeneities 
and effects of anisotropy, we have 
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(e. g., Rybicki & Lightmann 1979) where B := 5/(1 G). 

Evaluating the single-particle energy loss rate due to the 
SSC process, we restrict ourselves to regarding only single scat- 
tering events. If the particle distribution functions are given by 
power-laws, the synchrotron photon spectrum is also described 
by a power-law. However, since we are intrested in the detailed 
shape of the distributions, deviating from a simple power-law, 
we have to calculate the synchrotron spectrum in more detail. 

In an optically thin source, the differential number of syn- 
chrotron photons in the energy interval [e, e + de] is given by 

n sy (e) = 
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The summation symbol denotes the sum of the contribu- 
tions from electrons and positrons to the synchrotron spec- 
trum. Averaging over all pitch-angles of the electrons, the spec- 
tral emissivity P v (7) of a single electron of Lorentz factor 7 can 
be expressed as 
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(Crusius & Schlickeiser 1986) where a is the finestructure 
constant, 
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v e is the non-relativistic electron/positron gyrofrequency, 
v p is the plasma frequency of the relativistic pair plasma, i. e. 



7rm e (7) 



CS(x)=W *(x)W i(x) - Wib{x) Wis (a:) 



(21) 



(22) 



and W Ktll are the Whittaker functions. 

If the particle density exceeds n^lO 8 cm -3 synchrotron self 
absorption (SSA) is negligible for the following two reasons: 
The optical depth due to synchrotron self absorption for an 
initial power-law distribution 71(7) = no 7~ 2 can be estimated 
as 



TSSA 



y/ZctTlC 2 
m e c 2 



r(|)r (|) 



(23) 



(Schlickeiser & Crusius 1989). This optical depth becomes 
unity for 



r r, inW 1/3 D 2/3 1/3 r,l/3 -TJ 

^ssa = 5.7 • 10 7 1 / 3 B\ n g ' Hz 



(24) 



where 71 ; 3 is the lower cutoff Lorentz factor in units of 
10 3 , B-i is the magnetic field in units of 0.1 G, no is the pair 
density in units of 10 9 cm~ 3 and Rn is the blob radius in 
units of 10 11 cm. The frequency ^ssa is of the same order as 
the Razin-Tsytovich frequency 



^=2-10 11 n 9J BZ 1 ( 7 )- 1/2 s- 1 



(25) 



where the luminosity is strongly suppressed due to plasma 
effects. Here, (7) denotes the average electron Lorentz factor. 
This demonstrates that for pair densities ri>10 8 cm~ 3 syn- 
chrotron self absorption can be neglected. For lower densities 
we include it in our calculations, evaluating 



tssa = — Re 



00 

/ 



2 d ( n[f} 



d~f V 7 : 



(26) 



(e. g., Rybicki & Lightman, 1979) self-consistently. 

From the point of view of reacceleration, synchrotron self 
absorption can be neglected since particles resonating with 
photons at the lower cut-off of the synchrotron spectrum should 
have Lorentz factors of 



1/6 



7 res w 450 7 1 ' 3 ± J _ 1 



B' 



i/ei/e pi/6 



"9 



(27) 



which is lower than the particle Lorentz factors we deal 
with in the simulations carried out in this paper. 

If all scattering occurs in the Thomson limit, j± e <JC 1, the 
first order SSC energy loss rate is easily determined by 



(drft\ _ I 



car 



(7± - l) J deen sy (e). 




V dt J ssc 

where the synchrotron luminosity can be calculated as 

'd7± 



(28) 



J deen sy (e) 




Rb f 
3c J 



d7±n±(7i 



(drf±\ 
V dt J, 



(29) 



yielding for the initial power-law distribution functions 



(drn\ 
V dt / ssc 



7 



B 2 



16 2 

27 V 87rm e c ; 



r, 3-s s-1 s — 1 2 

2 , R B n± 7 2± 7 1± 7 . 

' o — s 



(30) G(E) 



Discarding the Thomson approximation and using the no- 
tation of the previous subsection, the exact expression for the 
SSC energy losses is 



to lowest order reduces Eq. (32) to Eq. (28) for j 2 > 1. In 
Fig. 4, the analytic solution (32) is compared to the Thomson 
scattering result (Eq. [27]). 



dj\ CTvr" 
ssc 2 

' 2 



dfj, / dn / de fn p (,(e) (1 — (3/j,) ■ 



1 +F(F- 1) + Fk u f 2 l-Pn + Pk' [/x - 



F 3 



- 1 



10' 
10 s 
10= 

T « 1 ° 4 

■&10 3 
10 2 
10' 
10° 



o = 4.0 










o = 4.5 


- 


" o = 5.0 





10" 



10 10 10 



Fig. 4. Energy loss rates of a test electron/positron due to in- 
verse-Compton scattering of synchrotron photons. B = 1 G, 
i?B = 10 15 cm, n e = 10 4 cm" 3 . Thomson-scattering (dotted curves) 
and exact Klein-Nishina cross-section (solid curves) 

In this case of an isotropic radiation field in the blob frame, 
the integrations over n and k can be solved analytically if we 
neglect terms of order I/7 2 <g; 1 and of order e/7 <§C 1. This 
yields 



00 

_ ( ~~Jt) =cirr 2 den ph (e) 
\at / ssc J 



G( 7 e) 



(32) 



where cur 2 G(^jt) /e cm -3 is the energy loss rate of an elec- 
tron of energy 7 scattering an isotropic, monochromatic radi- 
ation field of photon energy t and photon density 1 cm -3 and 



g( E ) = If 1 + 5E 



4:E 



3 (1 + 4E) 2 1+4-E 
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(33) 



For small values of E = ye (i. e. the Thomson limit) the 
Taylor expansion 
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Fig. 5. Energy loss rates of a test particle (Lorentz factor 74) 
due to inverse-Compton scattering of external (accretion disk) pho- 
tons (EIC; disk luminosity L = lO^ergs -1 , M BH = 1O 6 M , 
2 = 2- 10 -4 pc, T = 10) and of synchrotron photons (SSC; B = 1 
G, R B = 10 15 cm, n± 
emission (B = 1 G) 



10 4 cm 3 , a = 4.0), and to synchrotron 



It should be noted that for such an ultrarelativistic particle 
distribution Compton scattering its own synchrotron photons, 
Klein-Nishina corrections do not lead to a break in the energy 
dependence of the energy-loss rate, but due to the fact that 
with increasing energy a decreasing fraction of the scattering 
events occurs in the Thomson regime, lead to a much smoother 
flattening of the energy depencence than what is obtained 
from the relatively sharp soft photon distribution coming from 
the accretion disk. In the energy regime considered here, the 
energy-loss rate scales as —dy/dt ~ 7° with a w 1.5 for the 
spectral index s ~ 2 of the particle distributions n(j) ~ 7~ s - 
Clearly, this effect depends strongly on the temporal particle 
distribution which determines the synchrotron spectrum. 

The result of Eq. (32) (using Eq. [34] for E < 0.1) is in- 
cluded in Fig. 5 for a magnetic field strength of B = 1 G and 
a blob radius of Rb — 10 cm. Higher-order SSC scattering 
(up to nth order in the nth time step) is incorporated in our 
simulations by replacing n sy (e) by n sy (e) +nssc(e) in Eq. (32). 

3. Discussion of the elementary processes 

From Fig. 5 we can deduce some estimates about which of the 
various energy loss and exchange processes discussed above 
play an important role in different energy regimes of particles. 

First, we note that the energy transfer rates due to colli- 
sional effects (elastic scattering, bremsstrahlung emission) and 
the pair annihilation rate scale linearly with particle density 
n. The synchrotron and SSC losses basically scale as ~ B 2 
(the SSC losses, additionally, ~ Rb ■ n), whereas the exter- 
nal inverse-Compton losses basically scale as Lo / z 2 . The SSC 
losses also strongly depend on the steepness of the initial dis- 
tribution functions since a harder particle spectrum yields a 
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harder synchrotron spectrum which, in turn, leads to higher 4. The 7-ray photosphere 

inverse-Compton losses. 

The effect of elastic scattering is of minor importance for 
ultrarelativistic particles. We find the time-scale for thermal- 
ization to be given by 



Tth := 



It 



1 



!d 7t M| ( 



2.5 • 10 1 

ri4 



(35) 



for 74 = 100 which increases rapidly towards higher particle 
energies. Here, 714 is the particle density in units of 10 4 cm -3 . 

The distribution functions' high-energy tail will be most 
heavily influenced by one of the radiative energy-loss processes 
discussed above. Of course, the relative importance of the dif- 
ferent radiation processes is strongly dependent on the exact 
values and variation of parameters, especially the magnetic 
field strength and the luminosity of the central photon source. 

The time-scales for the respective processes can be esti- 
mated as 



TSSC := 



It - 1 

d<y t /dt\ 



5- 10 1 



ssc 



7 t 1/2 n 4 -Bo -Rii 



(36) 
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- where Rn - 
Compton process, 

3- 10 7 

for the synchrotron process, 



for the synchrotron-self- 
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(38) 



for the 



l() 46 crgs- 1 10-" pc 

external inverse-Compton losses in the Thomson regime, and 



_ 3.2- 10 11 _(,.! 

^"brcmsstr. ~ 'Tt S. 

724 



(39) 



In the beginning phase of the jet evolution, the effect of pair 
annihilation makes a negligible contribution to the distribution 
functions' evolution. The annihilation timescale for particles 
near the lower cutoff can be estimated as 



n(lt 



\dn{ lt )/dt\ 



10 

714 



12 



7t s 



(40) 



which indicates that pair annihilation and bremsstrahlung 
are of minor importance to ultrarelativistic pair plasmas of 
density ~ 10 9 cm -3 . 

These considerations and the investigation of the result- 
ing photon spectra lead us to some general simplifications in 
the treatment of the particle spectra evolution as long as the 
plasma is ultrarelativistic: 

(a) Pair annihilation is negligible. 

(b) Elastic (M0ller and Bhabha scattering) and inelastic scat- 
tering (bremsstrahlung emission) are only important for rela- 
tively low-energetic particles (7 < 100) which do not exist in 
our model system during the phase we are dealing with in this 
paper. 

(c) The effect of energy dispersion is neglected. 

(d) Bremsstrahlung emission is negligible for pair densities 



We now discuss the initial height of a plasma blob from which 
TeV 7-ray emission as observed from Mrk 421 can escape 
freely. Since BL Lac objects have very dilute emission line 
clouds, backscattering of accretion disk radiation into the tra- 
jectery of emitted 7 -rays is negligible. 
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Fig. 6. Pair production optical depth for high-energy 7-rays, emit- 
ted at different heights (zi), interacting with direct accretion disk 
radiation (L = 10 44 erg/s, M = 10 6 M ) 

Then, the optical depth for a photon of dimensionless en- 
ergy £1 , emitted at height zt due to interaction with accretion 
disk radiation is calculated as 

T 77 ((El,Zi) = 27T- 

1 OO 

dz J dn{l — n) J dea-,j(ei,e,fi)np h (e,fl;z) (41) 

Zi -1 2 

£1 (I-C) 

where n p h(e,f2; z) is the differential number of accretion 
disk photons at height 2, <r 77 is the pair production cross sec- 
tion 



cr 77 (/3) = 



16 



a T (l-p 2 ) [3-/3 4 ]ln 



1+13 



1-/3 



13 = 



eei (1 - /ti) 



2/3 [2-/3 2 ] ) , (42) 



(43) 



and we assumed for simplicity that the 7-ray is emitted 
along the symmetry axis. The result for different injection 
heights Zi is illustrated in Fig. 6. and shows that the 7-ray 
photosphere for 7-rays of energy 500 GeV - 1 TeV is located 
slightly above 2 • 10~ 4 pc in the case of an accretion disk of 
luminosity L — 10 44 erg s _1 . 

5. Simulations of the particle evolution 

Now that we know the energy-loss rates of the pairs inside a 
new jet component, we can follow the evolution of such a pair 
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plasma. In the first phase in which radiative losses are the only 
really important mechanisms, this can be done most easily by 
following the cooling of single particles and using the fact that 
the total particle number is conserved: Let 



f(l[t],t) = 



~dt 



be the cooling rate for one particle. From to to t 
the particle energy reduces from 70 to 74 where 

7t ~ 7o + At /(7o,*o)- 

Particle number conservation implies 
^70 



M7t) = n o{lo) 



d-yt 



(44) 
= to + At 



(45) 



(46) 



If we choose At small enough so that the explicit time de- 
pendence of / (which, e. g., represents the temporal variation of 
the background photon distributions) within one time interval 
is negligible, we find 



d'yo _ /(7o,M 
djt ~ f{lt,to) ' 



(47) 



This is used to follow the cooling of the injected pair 
plasma. In each time step, we calculate the emanating pho- 
ton spectra (see Sect. 6) and use the respective synchrotron 
and SSC spectra in order to calculate the SSC cooling-rate 
according to Eq. (32) self-consistently. 

6. Photon spectra 

From the observational point of view, it is most interesting to 
know both the evolution of the radiation spectra with time 
and the time-integrated photon spectra, which are emanating 
from the pair plasmas treated in the previous sections. We 
first evaluate all the contributions to the photon spectra in the 
frame of reference of the jet (blob frame) and then transform 
to the observer's frame using 



(48) 



N; h (e*,n*,t*) = D 2 N ph ( -,n,t 



where the asterisk denotes quantities measured in the 
observer's frame (which we identity with the accretion-disk 
frame), N(e,Cl) = Vb n{e, CI) is the spectral photon produc- 
tion rate in the energy interval [e, e + de] (e = hv / (m e c 2 )) in 
the solid angle interval [CI, CI + dCt], integrated over the whole 
blob volume Vb, D = (r[l — /3r/i*]) _1 is the Doppler factor, 
H* — cos 9* is the cosine of the observing angle in the observer's 
frame and /i = rr^A, is the observing angle cosine in the blob 



frame. A time interval A t measured in the blob frame is related 
to a light reception time interval A t r by A t r = A t / D and to 
an accretion disk frame time interval At* by At* = F At. 

The most important contributions to the 7-ray spectra in 
the MeV - TeV regime come from inverse-Compton interac- 
tions. The photon spectrum resulting from inverse-Compton 
scattering of accretion disk photons is calculated as 



h E ic^s,Cl 3 ) 



oc 



dCl c n ph (e, Cl e ) ■ 



d~f 



47T 



de s dCl s 



(49) 



where e s = hv a /(m e c 2 ) is the normalized energy and Cl s the 
solid angle of the motion of the scattered photons, n p h(e, Cl e ) 
is the differential number density of accretion disk photons at 
the location of the scattering event as discussed in Section 2.1 
and the angle variables are illustrated in Fig. 2 and Eqs. (6) 
and (7). 

It has turned out to be most convenient for this calculation 
(using the full Klein-Nishina cross section) to use the differen- 
tial Compton scattering cross section in the blob frame: 
d 2 a c 1 rl 



de s dCl s 



2 72 (1 - /3/x) 2 
7(1- fa) 



'e(l-K)+ 7 (!-&«.) 



e[l-/^] 



where 



3 [1 - P^s 



K =1- 



[i - PiA 



l+K 



/2 



(50) 



1-K 



= 1 + 



7 2 (1 - M (1 - fills) 
1 1 



7e(l-0A*) 7C.(1-/?A*.) 



(51) 



(Jauch & Rohrlich 1976). (Of course, using the KN cross sec- 
tion in the electron system and transforming the final photon 
state back to the blob frame leads to the same final expression 
for the spectrum of scattered photons. It is more straightfor- 
ward to use directly the cross section [50].) 

With the assumptions we made above, the SSC spectrum 
is isotropic in the blob frame, and 



nssc(e 3 ) 

47T 



de 



n sy (e) +n S sc(e) 

47T 



/ 



dCl, 



de s dCls 



(52) 



Technically, we include higher order SSC scattering by us- 
ing in each time step the SSC photon density nssc of the fore- 
going time step together with the synchrotron photon density 
as seed photon field for SSC scattering. 

Using the cross section (50), Eq. (52) can be rewritten as 



nssc(e s ) = ^ J d 7~^r- I < h ' I ''/'- / (l l>- 



n ph (Fe s ) yjl-Pfj) 
I- Pit. I N 



N 



7(1-/9 A*) 



-l+K 



(53) 



where p is the azimuthal angle between photon and elec- 
tron direction of motion around the direction of motion of the 
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scattered photon, the other angle variables are the same as in 
Section 2.1, k' is given by Eq. (52), n p h denotes the sum of 
synchrotron and previously produced SSC photons, and 



JV:=7(l-/3/i)-e.(l-/e), 
1-/3 Us 



F = 7- 



N 



+ v/i - K 2 yr 



fit cos p. 



(54) 
(55) 
(56) 



If pair bremsstrahlung can contribute significantly to the 
time-integrated photon spectra (which we find to be the 
case if r2>10 6 cm~ 3 ) we evaluate the rate of production of 
bremsstrahlung photons per unit volume, h br , in the comoving 
fluid frame as 



/ 



nbr(e)=4c J d7_n_(7_) / d7 + n+(7+)- 



:(e,7-,7+>7cm) 



where we use the differential cross section 



(57) 



(-) 



8 r a 7 C 
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. 2 'Jem ^cm 
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'Jem €c 



+ 



'Jem £c 



(58) 



(Alexanian 1968) in the center-of-momentum frame of the 
scattering particles, a is the hyperfinestructure constant, and 
the integration limits are given by 



e < 7±/3± 



and 



7c 



7± 



7± 



(59) 



The calculation of the synchrotron spectrum has been de- 
scribed in section 2.2. Pair annihilation radiation is negligible 
for an ultrarelativistic pair plasma. 

After having calculated the resulting synchrotron and 7-ray 
spectra, we have to check whether 7 -rays of energies ~ 1 TeV 
can escape the emitting region. For this purpose we calculate 
the optical depth due to 7-7 pair production of high-energy 7- 
rays interacting with the radiation produced in the blob. This 
calculation is carried out in the blob frame where we assume 
all contributions to be isotropic (which is not the case for the 
external IC component; but for a 7 -ray traveling with small 
angle with respect to the jet axis, the assumption of isotropy 
even overestimates the pair production optical depth) . The re- 
sults show that even in the first time step (which, of course, 
is the most critical one) the 7-7 optical depth does not exceed 
one for parameters suitable to fit the broadband spectrum of 
Mrk 421 (see section 8). Nevertheless, we include the small ef- 
fect of 7-7 absorption when calculating the emanating photon 
spectra which is the reason for the total emission (corrected for 
7-7 absorption) being slightly lower than the EIC component 
in Fig. 7 a. The injection of pairs due to 7-7 pair production 
(for details see Bottcher & Schlickeiser 1997) is negligible even 
in the case of relatively high density as assumed for Figs. 7 and 



7. General Results 

In this section we present some general results which we col- 
lected during a series of simulations with varying parameters. 

An interesting feature is that the time-averaged SSC emis- 
sion of a pair distribution with a lower cut-off 7i>,100 follows 
a power-law in the X-ray range with energy spectral index 
0.5<,a<0.7 which is consistent with the X-ray spectra of flat- 
spectrum radio quasars (FSRQs) and low-frequency peaked BL 
Lacs (LBLs) where the hard X-ray emission can be attributed 
to jet emission. The X-ray spectrum slightly hardens with in- 
creasing lower cut-off of the particle injection spectrum. 




E y [MeV] 



Fig. 8. Time-integrated broadband spectrum from the evolving blob 
illustrated in Fig. 7 (total 7-ray spectrum: solid; EIC: dashed; SSC: 
dot-dashed; bremsstrahlung: dotted; synchrotron: long dashed) 

Our simulations showed that assuming a particle spectrum 
extending up to 72 ~ 10 5 - 10 6 where some EIC scattering 
occurs in the Klein-Nishina regime, the EIC spectrum is harder 
than in the case of pure Thomson scattering where we recover 
the classical result n(e) ~ e ~( s + 2 )/ 2 (Dermer & Schlickeiser 
1993 a) above the break in the photon spectrum caused by 
incomplete Compton cooling. 

A harder time-integrated EIC spectrum (assuming the 
same spectral index of the injected pair distributions) can also 
result if cooling due to EIC is very inefficient (i. e. for a blob 
at large distance from the accretion disk) implying that the 
effect of dilution of the soft photon field according to its z~ 2 
dependence dominates the evolution of the EIC spectrum. 

If a second mechanism contributes to the cooling of 
the radiating pair population, the resulting external inverse- 
Compton photon spectra do not significantly differ from the 
case of a pure EIC model. Equally the assumption of a lower 
cutoff does not have an important effect on the high-energy 7- 
ray spectrum. The hard X-ray to soft 7-ray spectrum from EIC 
becomes harder with increasing lower cutoff, implying that in 
case of a high lower cutoff this part of the photon spectrum is 
always dominated by SSC. 

Time-averaged SSC-dominated 7-ray emission can not ac- 
count for strong spectral breaks Aa>0.5 in the MeV range as 
observed in several FSRQs, even if a high lower cut-off in the 
particle in the initial particle distributions is assumed. This 
fact is well illustrated by Figs. 8 and 10. In contrast, such a 
cut-off can produce sharp spectral breaks (Aa > 0.5) in the 
time-averaged EIC spectrum. 
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Fig. 7. Instantaneous broadband spectra from the blob of simulation 1. Parameters: 71 = 10 3 , 72 = 10 6 , s = 2.0, Rg = 5- 10 14 cm, n = 10 4 
cm- 3 , B = 0.05 G, L disk = 10 44 ergs" 1 , M BH = 10 6 M Zi = 5 • 10~ 4 pc, T = 20, (9 = 5°. Total emission (corrected for 7-7 absorption): 
solid; IC scattering of accretion disk radiation: dashed; SSC radiation: dot-dashed; synchrotron radiation: long dashed, (all quantities in 
the observer's frame) 



We find that even in the case of an accretion disk of very 
low luminosity and injection of the pairs occurring relatively 
high above the disk EIC scattering will give a significant con- 
tribution to the photon spectrum in the GeV - TeV range (see 
Fig. 10). If our assumption of a two-component 7-ray spectrum 
with neither SSC nor EIC being negligible is correct, we con- 
clude that the intrinsic GeV - TeV spectrum of blazars will in 
general not be a simple power-law extrapolation of the EGRET 
spectrum. 

The particle density plays a central role for the broad- 
band emission from ultrarelativistic plasma blobs. The 7-ray 
spectra of blazars can well be reproduced by dense blobs 
(n >• 10 4 cm~ 3 ). But in this case, escape of high-energy 7- 
rays without significant 7-7 absorption implies the presence of 
a very low magnetic field and a very weak synchrotron com- 
ponent. This problem is very severe for high-frequency peaked 
BL Lacs (as Mrk 421) where the 7-ray component can extend 
up to TeV energies at which they can be absorbed very effi- 
ciently by the synchrotron component if one assumes that both 
components are produced in the same volume. 

In the case of very high density (n> 10 7 cm~ 3 ) the cut-off 
to lower frequencies in the synchrotron spectra is not caused 
by synchrotron-self absorption but by the Razin-Tsytovich ef- 



fect. The same effect would also suppress synchrotron cooling 
below the Razin-Tsytovich frequency (Eq. [25]) which could 
consequently lead to a storage of the kinetic energy in pairs as 
long as the magnetic field does not change dramatically and 
the jet remains well collimated. As mentioned above, this ef- 
fect can only be important if the synchrotron component (also 
above the Razin-Tsytovich frequency) is very low compared to 
the 7-ray component during a 7-ray outburst. 



Then, a time lag between a 7-ray flare and the associated 
radio outburst can give information about the geometry of the 
jet structure. We suggest that this fact could be a major reason 
for the difference between FR I and FR II radio galaxies: If the 
jet is well-collimated over kpc scales, the Razin-Tsytovich effect 
prevents efficient cooling of the relativistic pairs to Lorentz 
factors lower than the Razin Lorentz factor, and the kinetic 
energy of the jet material can only be released in the radio 
lobes of FR II galaxies. If collimation is inefficient and the jet 
widens up rapidly after injection, its kinetic energy can be 
dissipated at short distances from the central engine, which 
would result in an FR I structure. 
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Fig. 9. Instantaneous 7-ray spectra from the blob of simulation 2. Parameters: 71 = 6 ■ 10 3 , 72 = 5 ■ 10 5 , s = 2.0, Rb = 8 ■ 10 15 cm, n = 0.2 
cm -3 , B = 0.07, B oc z~ 0A , L disk = 10 43 ergs -1 , M BH = 10 6 M , z t = 10 -2 pc, r = 20, 9 = 2°. Total emission: solid; IC scattering of 
accretion disk radiation: dashed; SSC radiation: dot-dashed; synchrotron radiation: long dashed, (all quantities in the observer's frame) 



8. Model fits to Mrk 421 

We can now use our code and the results summarized in 
the last section to reproduce the broadband spectrum of the 
first extragalactic object significantly detected in TeV 7-rays, 
Mrk 421. Mrk 421 is a BL Lac object at z = 0.031 which ex- 
hibited a prominent TeV flare in May 1994 (Kerrick et al., 
1995). In this flaring phase, a flux of (2.1 ± 0.3) ■ 10~ 10 pho- 
tons cm -2 s _1 above 250 GeV has been found. We assume 
Ho = 75km(s • Mpc) -1 , account for absorption by the inter- 
galactic infrared background radiation (IIBR) using the model 
fits by (Stecker & de Jager 1996) which do not differ very much 
from each other for photons of energies <1 TeV and neglect 
absorption by accretion disk radiation scattered back by sur- 
rounding clouds (Bottcher & Dermer 1995). 

Our first example, shown in Figs. 7 and 8, demonstrates 
that jet parameters are possible where the 7-ray emission from 
Mrk 421 in its quiescent state can emerge from a small region 
where light-travel time arguments allow for even more rapid 
flickering in 7-rays than observed recently by Gaidos et al. 
(1996). As illustrated in Fig. 11, the synchrotron emission from 
this emission region is far below the observed radio to optical 
flux from Mrk 421. 

A successful fit to the broadband spectrum of Mrk 421 in 
its quiescent state can be achieved with a very low density, as 



chosen for our second example, illustrated in Figs. 9 and 10. 
This is shown in Fig. 12. 

The flaring state of Mrk 421 can result from an increase 
of the maximum Lorentz factors of the pairs which could be 
related to some increase of energy input in hydromagnetic tur- 
bulences accelerating the primary particles which can plausibly 
also imply an increase in the magnetic field. Increasing the cut- 
off energies to 71 = 8 • 10 3 , 72 = 7 ■ 10 5 and the magnetic field 
to B — 0.08 G while the other parameters remain the same is 
in our model calculation for the low state (see caption of Fig. 
9) yields an acceptable fit to the flaring state of Mrk 421, as 
shown in Fig. 13. Here, we integrated over an observing time 
of ~ 2 ■ 10 3 s, after which the the flux from this blob is far 
below the quiescent one. 

Petry et al. (1996) found a spectral index of the differential 
spectrum of Mrk 421 above 1 TeV of 3.6 ± 1.0 which translates 
to a spectral index of the integral spectrum of 2.6 ± 1.0. In 
Figure 14 we demonstrate that our model fit to the low state 
of Mrk 421 is in agreement with this measurement. Above 1 
TeV, absorption by the IIBR becomes important. In order to 
account for this, we used the model fits by Stecker & de Jager 
(1996). The figure demonstrates that both fits result in very 
similar absorbed spectra. The HEGRA flux value is obtained 
in averaging over long observing times compared to the vari- 
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Fig. 10. Time-integrated 7-ray spectrum from the blob illustrated 
in Fig. 9 (total 7-ray spectrum: solid; external IC: dashed; SSC: 
dot-dashed; synchrotron: long dashed) 
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Fig. 11. Comparison of the time-averaged photon flux from sim- 
ulation 1 (Figs. 7, 8) to the broadband spectrum of Mrk 421 in 
its quiescent state, indicated by squares (data from Macomb et al. 
[1996]) 
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Fig. 12. Comparison of the time-averaged photon flux from simu- 
lation 2 (Figs. 9, 10) to the broadband spectrum of Mrk 421 in its 
quiescent state 
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Fig. 13. Model fit to the broadband spectrum of Mrk 421 in its 
high state (stars); 71 = 8 ■ 10 3 , 72 = 5 ■ 10 5 , B = 0.08 G; other 
parameters as given in the caption of Fig. 9 

ability of Mrk 421. It lies between the Whipple data points for 
quiescent and flaring state (included in Figs. 11 - 13) and is 
well in agreement with a duty cycle of 5 % between high and 
low TeV state (Petry, priv. comm.). 
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Fig. 14. Integral photon flux resulting from simulation 2. The dotted 
and dot-dashed lines represent a power-law integral spectrum with 
spectral index 2.6 and 3.6, respectively. The solid and dashed lines 
are calculated describing absorption by the intergalactic infrared 
background radiation by Stecker & de Jager's (1996) model fit 1 
and 2, respectively 

It is remarkable that the broadband fits shown here imply 
an injection height of z; m 10 -2 pc which is far outside the 
7-ray photosphere for TeV photons. A high low-energy cutoff 
in the electron spectrum is required in order not to overpro- 
duce the X-ray flux. This, in turn, implies a very low particle 
density in order to recover the observed relative luminosity of 
the synchrotron to the SSC component. 

9. Summary and conclusions 

We gave a detailed discussion of the various radiative energy 
loss mechanisms acting in relativistic pair plasmas and concen- 
trated on the application to jets from active galactic nuclei. In 
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the first phase after the acceleration or injection of pairs radia- 
tive cooling due to inverse-Compton scattering is the dominant 
process. Here, Klein-Nishina effects are much more pronounced 
for inverse-Compton scattering of external radiation (from the 
central accretion disk) than for the SSC process since the syn- 
chrotron spectrum extends over a much broader energy distri- 
bution, allowing for efficient Thomson scattering at all particle 
energies. 

We developed a computer code which allows to follow the 
evolution of the energy distributions of electrons/positrons self- 
consistently, accounting for all Klein-Nishina effects and for the 
time- dependence of the synchrotron, SSC and external radia- 
tion fields which contribute to the radiative cooling of the pairs 
and computing the instantaneous 7-ray production from the re- 
sulting pair distributions. Since the intrinsic cooling timescales 
are much shorter than the time resultion of present-day instru- 
ments, we calculated time-integrated spectra and compared 
them to observations. 

We found some general results from a series of different 
simulations: 

(a) Due to the short cooling times involved in the SSC ra- 
diation mechanism, implying that we only see time-integrated 
SSC emission leading to an overprediction of hard X-ray flux 
from blazars, we favor EIC to be the dominant radiation pro- 
cess in 7-ray blazars. 

(b) Detailed simulations show that the EIC spectrum of 
a cooling ultrarclativistic pair distribution with maximum en- 
ergies implying EIC scattering in the Klein-Nishina regime is 
harder in the case of pure Thomson scattering where we recover 
the classical results of Dermer & Schlickeiser (1993 a). 

(c) The low-energy cutoff in the synchrotron spectra 
from dense pair jets is not determined by synchrotron self- 
absorption, but by the Razin-Tsytovich effect. Only when the 
jet widens up, radio emission at is<W n Hz can escape the jet. 
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